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The ion-bombardment induced evolution of intermixing is studied by molecular dynamics simu- 
lations and by Auger electron spectroscopy depth profiling analysis (AESD) in Cu/Co multilayer. 
It has been shown that from AESD we can derive the low-energy mixing rate and which can be 
compared with the simulated values obtained by molecular dynamics (MD) simulations. The overall 
agreement is reasonably good hence MD can hopefuly be used to estimate the rate of intermixing 
in various interface systems. 



Low-energy ion beams are commonly used in surface 
analysis and for film growth [1,2]. The use of ion- 
sputtering in the controllable production of nanostruc- 
tures and self-assembled nanoppaterns have also become 
one of the most important fields in materials science 
[2-5]. 

As the dimensions of the surface features is reduced 
to the nanoscale, many classical macroscopic (continuum 
and mesoscopic) models for morphological evolution lose 
their validity. Therefore the understanding of the driv- 
ing forces and laws governing mass transport involved 
in the synthesis and organisation of nanoscale features 
in solid state materials is inevitable. Unfortunately the 
fundamental understanding and the nanoscale control 
of atomic transport processes, such as interdiffusion is 
not available yet [6,7]. In order to get more insights in 
the atomic relocation processes during postgrowth low- 
energy ion-sputtering, it is important to measure and to 
calculate accurately the rate of intermixing at the inter- 
face. 

In this Letter we will show that using a newly devel- 
oped code for simulating ion-sputtering based on molecu- 
lar dynamics we are able to get mixing rates. We will also 
show that it is possible to extract the mixing rate data 
from AESD. It turns out that the agreement of the two 
methods is reasonably good. In this way we also could 
check the reliability of molecular dynamics simulations. 

The AESD [8,9] measurements were carried out on a 
Cu/Co multilayer system. The sample was made by sput- 
ter deposition on polished single-crystal (111) silicon sub- 
strates in a plasma beam sputter deposition system. It 
was characterized by XTEM and RBS and flat interfaces 
have been found [9]. AESD depth profiling was carried 
out using a dedicated device [10] by applying Ar + ions of 
energy of 1 keV and angle of incidence of 10° with respect 
to the surface of the crystal (grazing angle of incidence). 
The sample was rotated during ion bombardment. 

The depth profiles were measured as a function of the 
sputtering time keeping the bombarding ion current con- 
stant. A STAIB DESA 100 pre-retarded CM A with fixed 
energy resolution was used to record the AES spectra. 
The following AES peaks were detected Cu (60 eV), Cu 
(920 eV) and Co (656 eV, to avoid overlapping). A part 



of a typical depth profile is shown in Fig 1. For clarity 
only the copper is shown in Fig 1. To calculate the con- 
centration the intensity of the copper Auger current was 
normalized to that of the pure copper. The depth scale 
was calculated from the known thickness of the sample 
[9]. 

It is generally supposed that during ion mixing the 
atomic movements are similar to that of the usual diffu- 
sion and thus the occurring broadening can be described 
by the same equations [11]. Accordingly in case of a bi- 
layer system the concentration distribution formed due 
to ion bombardment of the interface by a given fluence 
can be described by the erf function. The variance, a 2 , 
of the erf function determines the extent of ion mixing. 
In many cases a 2 depends linearily on the ion fluence 
<f>. In these cases the <r 2 /<I> ratio called as mixing rate 
characterizes the mixing for a given ion bombarding con- 
dition. Important advantage of using the term mixing 
rate is that it can be in principle directly derived from 
the experiments. 

In the case of AESD we measure the Auger current 
of the elements present. The measured Auger current 
of element j Ij depends on the in-depth distribution as 
follows: Ij = + ij(2)nj + ij(3)n 2 + where ij(k) 
is the Auger current of element j emitted by the atomic 
layer k, and Kj gives the attenuation of the Auger current 
crossing an atomic plane. ij(k) ~ Xj(k)N(k), where 
Xj(k) is the concentration of element j in layer k, while 
N(k) is the number of atoms of the fc-th atomic plane. 
It is evident that in general from a single measured Ij 
one cannot determine the ij(k) values. On the other 
hand this equation can be used to simulate the measured 
Auger current during the depth profiling procedure if we 
assume an in-depth distribution. 

We do not know any experimental measurement of the 
in-depth distributions formed during AESD applying low 
(0.2 — 2 keV) ion energy. It seems, however, that dy- 
namic TRIM simulation can reliably be used to describe 
AESD [12]. This calculation also provides the in-depth 
distribution during the procedure. It turns out that at 
the beginning of the depth profiling procedure, when the 
interface is still far from the surface the in-depth distri- 
bution can be approximated by the erf function. For the 
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FIG. 1. The measured depth profile obtained at 1 keV ion 
energy. The concentration of Cu is given as a function of 
the removed layer thickness (measured from the interface in 
A). Empty diamonds denote the points at which a 2 values are 
determined by AESD. Inset figure: The broadening (a 2 ) of the 
interface as a function of the fluence (ion/ A ). The dotted 
lines show the two extreme of the possible slopes (mixing 
rates) . 




number of ions 

FIG. 2. The simulated filtered (R 2 ) (R2) as a function 
of the number of the ion impacts (number fluence) at 1 
keV ion energy at grazing angle of incidence using simulated 
ion-sputtering (molecular dynamics). Dotted line is a linear 
fit to the {R 2 ) curves. (R 2 ) is simulated with ion impact depth 
positions at the free surface (9 ML above the interface) and 
also at 4 ML above the interface. Inset figure: The simulated 
filtered (R 2 ) for sputter removal up to the removal of 4 layers 
as a function of the ion dose (ion/A 2 ). The error bars are 
standard deviations. 



evaluation of the experimentally determined depth pro- 
file we will suppose that (i) the interface has an intrinsic 
surface roughness, which can also be approximated by 
an erf function of er - (ii) ion mixing is the only process 
contributing to the broadening of the interface at least in 
the beginning part of the depth profiling, which will be 
studied. Thus for any measured copper Auger current we 
should find x (the distance of the interface from the sur- 
face) and a m (the measured variance) values of the erf 
function. Then the variance due to the ion bombardment 
induced mixing is <r 2 = — a\ . We derived a 2 at depths 
as indicated in Fig 1. at 1 keV ion energy at ~ 10° angle 
of incidence with respect to the surface. To proceed we 
must know the number of ions causing the broadening at 
the interface. In our experimental arrangement we can- 
not measure the ion fluence. On the other hand we can 
measure accurately the removed layer thickness. Taking 
the sputtering yield from the literature Y « 1.2 [9], we 
can derive the curve a 2 vs. fluence $ which is shown 
in inset Fig. 1. Inset Fig. 1 also shows the two limiting 
slopes; thus we can derive the mixing rate being 250 ± 150 

A 4 . 

Classical constant volume molecular dynamics simula- 
tions were also used to simulate the ion-solid interaction 
(ion-sputtering) using the PARC AS code [13]. Further 
details are given in ref. [13,15,16]. We irradiate the bi- 
layer Cu/Co (9 monolayers, (ML) film/substrate) with 
1 keV Ar + ions repeatedly with a time interval of 5-20 
ps between each of the ion-impacts at 300 K which we 
find sufficiently long time for the termination of interdif- 
fusion, such as sputtering induced intermixing (ion-beam 
mixing) [16]. The initial velocity direction of the imping- 
ing ions were 10° with respect to the surface of the crystal 
(grazing angle of incidence) . 



To describe homo- and heteronuclear interaction of Cu 
and Co, the Levanov's tight-binding potentials are used 
[17]. The cutoff radius r c is taken as the second neighbor 
distance. 

We randomly varied the impact position and the az- 
imuth angle </>. In order to approach the real sputter- 
ing limit a large number of ion irradiation are employed 
using simulations conducted subsequently together with 
analyzing the history files (movie files) in each irradiation 
steps. In this article we present results up to 100 ion ir- 
radiation which we find suitable for comparing with low 
to medium fluence experiments. 100 ions are randomly 
distributed over a 40 x 40 A 2 area. The size of the sim- 
ulation cell is 100 x 100 x 75 A 3 including 62000 atoms. 
The computer animations can be seen in our web page 
[14]- 

The simulation of the sputter-removal has been carried 
out in the following way: First we ion-bombarded the 
first (top) layer (the distance from the interface is -9 ML) 
using 100 ions, than using the ion-sputtered simulation 
cell again ion-sputtering is simulated by 100 ions in the 
2nd layer (-8 ML). Following this procedure in the 3rd, 
4th and 5th layers we get an ion-sputtered sample with an 
intermixed interface. We find that (R 2 ) remains rather 
small until the ions are initilaized in the 3-4th layers (- 
5 ML, see Fig. 2). This is because the distance from 
the interface is beyond the projected range of the 1 keV 
ions in the upper layers. To reduce the computational 
demand, we do not remove completely the layers. 

In Fig 2 the evolution of the sum of the square atomic 
displacements (R2) along the depth direction of all in- 
termixing atoms (R 2 ) = ^2f[zi(t) — Zi(t = 0)] 2 , where 
Zi(t) is the depth position of atom i at time t, can be fol- 
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FIG. 3. The simulated broadening (in ML) as a function of 
the distance from the interface (in ML) at 1 keV ion energy at 
10° angle of incidence obtained by simulated sputter-removal. 
—9 ML corresponds to the free surface. The error bars are 
standard deviations. 



lowed as a function of the number of ions. Using the rela- 
tion a 2 = (R 2 ) it is possible to calculate the mixing rate 
k = (R 2 ) I $ from the slope of the fitted straigh line on the 
(R 2 ) vs. $ curve in Fig 2. The contribution of intermixed 
atoms to (R 2 ) is excluded (filtered (R 2 )) in a given layer 
if the in-layer concentration of them is less then ~ 16 
%. The filtered simulated (R 2 ) is used to compare with 
the measured a values. The lower concentratioan limit 
of ~ 16 % corresponds to the usual 84 — 16 rule used for 
extracting the depth resolution in AESD [19]. Therefore, 
the comparison with AESD is meaningful when contribu- 
tions to broadening (and to (R 2 )) are filtered out which 
gives less in-layer concentartion than ~ 16 %. The re- 
sults of simulated sputter-removal is shown in the inset 
Fig 2. We also show results in Fig 2 for simplified simula- 
tions, when ions are directly initialized from between the 
4th and 5th layers (embedded ion-sputtering) using also 
grazing angle of incidence. The evolution of the filtered 
(R 2 ) as a function of the number of ions shows us that 
(R 2 ) is very close to that obtained for simulated sputter 
removal. Therefore we conclude from these results that it 
is sufficient to simulate intermixing with embedded ion- 
sputtering. 

The calculated mixing rate on the basis of inset Fig 2, 
is k w 350 ± 100 A . This value is reasonably compara- 
ble with the value obtained by AESD measurements. It 

must also be noted that Cai et al. obtained k « 400 A 
using 1 MeV Si + ions and X-ray scattering techniques 
in Co/Cu multilayer [18]. Since the ion mixing occur at 
the low energy end of the cascade process, e.g. the high 
energy collisional cascades split into low energy subcas- 
cades (see e.g. refs. in [15]), an agreement is expected 
for. We also give the measured and calculated mixing ef- 
ficiencies of £ = k/Fo [11,15], where Frj is the deposited 
ion-encrg y/depth (F D « 45 cV/A obtained by SRIM [20] 

simulation) . We get £ « 6 ± 3 A / eV by experiment and 

£ « 7 ± 2 A /eV by simulations. On the basis of this 
value we can characterize low-energy ion-sputtering in- 



duced intermixing in Cu/Co as a ballistic interdiffusion 
process. 

The saturation of broadening is reached at around ~ 
9 ± 1 ML interface width as it is shown in Fig 3 obtained 
by simulated sputter-removal. 

In this Letter we have presented that the combina- 
tion of atomistic simulations with Auger electron spec- 
troscopy depth profiling might be a new efficient method 
to depth profiling analysis of multilayered materials. 
Also, the reasonably good agreement between experiment 
and simulations provides us the possibility of predicting 
interdiffusion properties for various multilayers for which 
no experimental results are available. 
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